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Distributions of Coulomb blockade peak spacing are reported for 
large ensembles of both unbroken (magnetic field B = 0) and broken (B ^ 
0) time reversal symmetry in GaAs quantum dots. Both distributions are 
symmetric and roughly gaussian with a width ~ 2-6% of the average 
spacing, with broad, non-gaussian tails. The distribution is systematically 
wider at B = by a factor of ~ 1.2 ± 0.1. No even-odd spacing correlations 
or bimodal structure in the spacing distribution is found, suggesting an 
absence of spin-degeneracy. There is no observed correlation between peak 
spacing and peak height. 
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For some time it has been appreciated that electron transport in mesoscopic systems 
exhibits quantum interference effects with universal statistical features, and that this 
universality can be associated with the underlying universality of quantum chaos [1] and its 
mathematical description in terms of random matrix theory (RMT) [2]. This approach has 
been successful in describing transport in open quantum systems (i.e. systems with large 
conductance, g>e^ Ih, to reservoirs) where a single-particle picture apparently provides 
an adequate description of the physics. Recent application of RMT to ground state 
properties of nearly isolated quantum dots — in particular, in characterizing the distributions 
of Coulomb blockade (CB) conductance peak heights [3] has also been remarkably 
successful [4]. 

On the other hand, experiments by Sivan et al. [5] and Simmel et al. [6] suggest 
that the most basic prediction of RMT, namely the famous Wigner surmise for the 
distribution of level spacings, fails to describe the fluctuations of Coulomb blockade peak 
spacing, implying that fluctuations in the energy separation between adjacent ground states 
of a quantum dot — the so-called addition spectrum — appear not to be distributed according 
to RMT. In particular, these experiments [5, 6] found CB peak spacing fluctuations of 
order 0.1-0.15 of the average spacing, larger than predicted by RMT assuming constant 

charging energy E^ = e^ I C^^^ with the total capacitance of the dot given by C^^, . The large 
fluctuations observed in the experiment and supporting numerics lead Sivan et al. [5] to 
suggest that classical charging energy fluctuations proportional to E^. not included in RMT 
dominate peak spacing fluctuations. On the other hand, recent random phase 
approximation (RPA) calculations lead to the opposite conclusion, that fluctuations due to 
charge rearrangement should be smaller than [7] or of order [8] the mean single-particle 

level spacing A, much smaller than E^,. However, one should use caution in applying 
these theories to semiconductor quantum dots, since RPA breaks down at low electron 
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densities where single-particle kinetic energies and interparticle potential energies are 
comparable. 

In this Letter, we present an extensive experimental study of the spacings of CB 
peaks in GaAs quantum dots, for both zero and non-zero magnetic field, including over 
20,000 CB peaks measured in seven devices. We find that the distributions of CB peak 
spacing fluctuations are not well described by single-particle, spin-degenerate RMT 
(SDRMT); the main discrepancy being the absence of bimodal structure in the measured 
distributions. Both B = and B ?t peak spacing distributions are roughly gaussian with 
non-gaussian tails, in qualitative agreement with [5, 6]. In contrast to previous 
experiments, however, we find that the width of the peak spacing distribution is narrow, 
comparable (once scaled) to the mean level spacing, and shows the effects of time reversal 
symmetry breaking comparable to RMT predictions, suggesting that quantum effects play a 
role in determining the distributions. We also find no correlation between CB peak heights 
(reflecting eigenfunction properties) and peak spacings (reflecting eigenvalue properties). 

Ground state energy fluctuations are measured using the Coulomb blockade, which 

appears in quantum dots with tunneling leads (left, right lead conductance gi,gr< 2e^ Ih), 
when temperature T and source-drain voltage V^^ are less than E^. [9, 10]. In this regime, 
dot conductance is suppressed except when the (A'-i-l) and A'^ electron ground state energies 
of the dot differ by the chemical potential of the leads; this degeneracy condition, when 
tuned by a gate voltage , produces a series of narrow conductance peaks nearly periodic 

in . At very low temperature and bias, ( k^T, eV^^ ~< A (we define A = iTZYi Im A as 

the average single particle level spacing for spin-degenerate levels; m* is the electron 
effective mass, and A is the area occupied by the electrons), transport on a CB peak is a 
resonant process, making the peak position sensitive to the discrete spectrum of the dot. A 
simple model connecting level spacing statistics and CB peak spacing is the so-called 
'constant interaction' (CI) model, in which the separation between ground state energies of 
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the dot probed by CB is separated into two parts, a charging energy E^. that is independent 
of (or at most slowly varying with) the number of electrons on the dot, A^, and a fluctuating 
part associated with a discrete quantum level spacing. This separation assumes that 
fluctuations in due to charge rearrangement in the ground state upon adding an electron 

are small compared to A. Within the CI model, the spacing (in gate voltage) between CB 
peaks, AV' = V'^^ - V' where V' is the center position of the peak, is given by 



where e' is the single-particle energy level and 7]= C^l C^^f is the ratio of gate 

capacitance to total dot capacitance, also assumed to vary slowly with /. The dependence of 
Eq. 1 on whether i is even or odd reflects the spin degeneracy of levels. 

If one further assumes that spacings (e'^^ -e') obey spin-degenerate (SD) RMT 
statistics, the resulting distribution P{v) of normalized fluctuations in peak spacing v = 

^AV^' -^Ay^'^jy^^AV^'^ will consist of a 5 function for the case of i odd plus a Wigner- 

Dyson distribution (P^^ (5) oc se'""'^* for B = and oc e"^"'/" for B ^ 0, with s in units 
of A) for the case i even. The brackets (•) denote an average over an ensemble of peaks. 
Note that within this CI+SDRMT model, the standard deviation, tr , of P{v) depends on 
time-reversal symmetry; for A « E^, CI+SDRMT gives G =0.62(0.58) A /E^ for B = 
(B ^ 0) and a ratio Gg^^ I g^^q ~ 1 . 1 . For GaAs dots of the type described here and in Ref . 
[5, 6], the predicted peak spacing distribution width is g{v) ~ 0.03, depending slightly on 
the dot shape and size, based on experimental values of A and E^ (see Table 1). 




(z even) 
(/ odd) 



(1) 
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Since a bimodal P{v) is not found experimentally, as seen below, one is motivated 
to consider a CI model without even-odd structure, based on spin-resolved (SR) level 
statistics. The breaking of spin degeneracy, discussed in [7] and observed experimentally 
in [11], will lead to level statistics described roughly by two overlapping Wigner-Dyson 
distributions [12]. The ratio (7^^0/(75^0 = 0.70/0.65-1.1 in CI+SRRMT is similar to 

CI+SDRMT; however, the predicted spacing distribution width (7(v) ~ 0.01 is reduced by 
a factor of ~ 2 due to a reduced mean level spacing, A^^j = A/2 . 

The quantum dots were formed by gate depletion of a two-dimensional electron gas 
(2DEG) in a GaAs/AlGaAs heterostructure (see Table 1). Irregular dot shapes were 
designed to produce chaotic scattering, and all devices were smaller than the bulk mean free 
path so that transport within the dots is ballistic. The dot area. A, was estimated from the 
lithographic area with a depletion of -120 nm around the perimeter of the device. Charging 

energy was measured from the relation E^. ~ erjl^AVgj, based on Eq. 1 with A « E^, with 

77 extracted from the T dependence of the peak width [10]. AH measurements were made 

using two- wire ac lockin techniques with a voltage bias of 5 |i.V at 11 Hz. The electron 

base temperature, determined by fitting peak width versus temperature [10], was -100 mK 
for all devices. Ensemble statistics were collected by sweeping one gate voltage, , over 

- 20 peaks then incrementing magnetic field or a second gate voltage to yield a new 
ensemble of peaks. All data are based on relatively small CB peaks, in the range 0.01 to 

0.1 e^/h. 

A t5^ical scan of CB peaks is shown in Fig. la. To extract peak spacing, each 
peak is fit by a cosh'^ form [10] and spacing is determined from the centers of the fits (Fig. 
lb). Since the average spacing decreases with increasing A'^ (reflecting an increasing C^^,), 

a running average (AV^^ is found from the best fit line (dashed in Fig. lb) and used to 
define V = (aV^' -^AV^'^jy^^AV^'^. Experimental noise in the spacing distribution (for 
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instance due to charge rearrangement in the doping layer), given as <Tf^„i,^(v) in Table 1, 
can be separated from real spacing fluctuations by comparing measurements at opposite 
magnetic fields (Fig. Ic). The noise c^voj^gCv) is defined as the distribution width of 
fluctuations antisymmetric in B and so would not include field independent gate voltage 
induced charge rearrangement. While the heights of nearby CB peaks show considerable 
correlation (Fig. la), peak spacings appear uncorrelated (Fig. lb, Ic). We estimate the 
number of "independent" peak spacings, n. in Table 1, as the number of peaks measured in 
each scan of multiplied by the number of peak scans with characteristically different 
heights. 

Histograms of peak spacings from three dots with similar device parameters (dots 
3, 4, and 5) are shown in Fig. 2. Both the B = and B =^ histograms are roughly 
symmetric and gaussian. A gaussian fit to the B = histogram gives a standard deviation 
0'b=o(^)/ = 0.019 (the subscript/ indicates 'gaussian fit') whereas a direct evaluation of the 
second moment of the spacing data set yields ^^^^(v) = 0.027. This difference results 
from broad non-gaussian tails, which can be seen a logarithmic plot (right insets. Fig. 2). 
All B data are taken with between 3-15^0 through the device, where 
% = h/e = 4.14 mT[Lm^. The B ^ distribution width from a gaussian fit yields 
<Tg^Q(v)f = 0.015 and the second moment of the data gives <Tjj^q(v) = 0.022. By either 
method, a ratio Ob=qIOb^q ~ 1.2-1.3 is obtained. From aU devices, 
(Tg^o/(Tj,^o = 1.2±0.1 (see Table 1), slightly larger than but comparable to the value 
predicted by either CI+SDRMT or CI+SRRMT. The left insets in Fig. 2 show that the 
CI+SDRMT spacing distributions for these devices convolved with (gaussian) 
experimental noise remain bimodal, clearly inconsistent with the data. 

The B^Q peak spacing histogram for the smallest, quietest device (#1) is shown in 
Fig. 3a along with a no-adjustable-parameters CI-i-SRRMT distribution convolved with 
experimental noise. An important consideration is how the width of this distribution is 
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affected by temperature. Figure 3b shows <yB*o{^B^/^sR) measured for two gate voltage 
configurations of device 1 (with different A^^'s) along with a simple model based on 
thermal averaging of the peak spacings. In the model, finite-temperature peak spacings 
= 'LjWi j Sj are constructed from zero-temperature spacings £') obeying 

SRRMT statistics with (5,) = A^^ . The weighting function w^ j <x cosh-^{ \e' - /2kjJ') is 
chosen to include more contributions from adjacent spacings 5, at higher temperature (not 
including height fluctuations) leading to <t(T) ocT-'" in the limit fc^T » A^^, consistent 
with the data. While this model appears to have the correct cross-over scale and correct 
limits at fcgT « A^.^^ and k^J^ » A^j^, a proper theory of cr(T) has yet to be developed. 
However, the observed temperature dependence suggests a quantum rather than classical 
origin to the spacing fluctuations. 

Finally, we have looked for correlations between peak spacing fluctuations and 
peak height fluctuations. Within RMT, no correlations are expected since height 
fluctuations depend on eigenfunctions while spacing fluctuations depend on eigenvalues 
which are uncorrected with eigenfunctions [13]. Experimentally, we find that the 

normalized spacing v between adjacent peaks (i+l and i) is not correlated with either the 
normalized height difference, Ag = (g'^l^- g'^^)/{^g^^), or their normalized average 

height, g = (g'^.l + gLj/2(g„,,), as shown in Fig. 3c and 3d. 
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Captions 



Fig. 1. (a) Coulomb blockade peaks (diamonds) at B = 30 mT as a function of gate voltage 
for device 1. Solid curve shows fits to cosh'^ lineshape. Left inset: Detailed view of 

data and fit on log-linear scale. Right inset: Micrograph of device 1 , other devices are 
similar, (b) Peak spacings extracted from data in (a) at B = +30 mT (diamonds) B = -30 

mT (open circles). Dashed line is best fit (to +30 mT data), corresponding to (AV^'^. (c) 

Dimensionless peak spacing fluctuations, v = {^AV^ - ^AV^'^j/^AV^'^, as a function of gate 
voltage Vg for data in (b). Differences between ±30 mT data indicated experimental noise. 
Normalized (spin-resolved) mean level spacing A^j^l E^, indicated by vertical bar (See 
Table 1). 

Fig. 2. Histograms of normaUzed peak spacing v (bars) for (a) B = and (b) B for 

devices 3, 4, and 5. Solid curves show best fit to normalized gaussian of width 0.019 
(0.015) for B = (B ^ 0). The B = histogram is wider by a factor of ~ 1.2 than the B ^ 
histogram. Data represents 4,300 (10,800) CB peaks from the devices with ~ 720 
(1600) statistically independent for B = (B 0). Horizontal bar indicates (spin-resolved) 
mean level spacing A^^j / averaged over the 3 devices. Right insets: Plots of histogram 
(diamonds) and best fit gaussian (solid curve) on log-linear scale. Dashed curve is 
gaussian of width 0.13 from Ref. [5]. Left insets: Dotted curves are CI+SDRMT peak 
spacing distributions; solid curves correspond to CI+SDRMT distributions convolved with 
gaussian of width C7^„,jg(v) = 0.009 averaged over the three dots (see Table 1). 

Fig. 3. (a) Histogram of peak spacing data from dot 1 at B (bars). Dotted curve 
shows zero-temperature B ^^^^ CI+SRRMT peak spacing distribution; solid curve is 
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CI+SRRMT distribution convolved with gaussian of width (T^^,.^^(v) = 0.006. (b) 
Temperature dependence of cjg^Q(v) for two gate configurations of dot 1 ( A^^ = 21 |j,eV 
solid diamonds; A^,^ = 14 |i.eV open circles), and thermal averaging model (solid curve) 
described in text, with A^jj = 21|j,eV and E^ = 590[LqV. (c),(d) Scatter plot of 

normalized peak spacing v versus adjacent peak height differences 

= (g'r^ - gL)/(gmax) and avcragcs g = [g'^ + gUj/2{g^) for all B data from 
dot 4. No correlations between spacings and average heights or height differences are 
observed. 
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Parameter 


dot 1 


dot 2 


dot 3 


dot 4 


dot 5 


dot 6 


dot 7 


A (|LLm^) 


0.17 


0.20 


0.32 


0.34 


0.38 


0.47 


0.50 


o 

d(A) 


900 


900 


800 


800 


900 


800 


900 


(^eV) 


21 


18 


11 


11 


9 


8 


7 


Ec (l^eV) 


590 


760 


580 


500 


380 


600 


320 




340 


400 


900 


1000 


800 


1400 


1000 


n 

I 


190 


70 


140 


830 


1300 


710 


420 


B=0 




48 


38 


25 


25 


43 


56 


(xlO~^) 




(V) 


(2) 


(2) 


(3) 


(2) 


(3) 




18 


34 


23 


22 


20 


38 


43 


(xlO~3) 


(2) 


(4) 


(3) 


(1) 


(2) 


(2) 


(2) 






1.3 


1.7 


1.2 


1.2 


1.1 


1.3 






(0.2) 


(0.2) 


(0.1) 


(0.2) 


(0.1) 


(0.1) 


Noise^ ^ 


6 


23 


15 


10 


8 


25 


30 


(xlO~3) 

















Table 1. Device parameters and measured spacing statistics at T ~ 100 mK: area occupied 
by electrons (A), 2DEG depth (d), SR mean level spacing (A^^ = Ttli' / m*A), charging 

energy {E^-e^ I number of electrons in dot (AO, number of statistically independent 
peak spacings {n), and peak spacing distribution width at B = (c^^qCv)) \^ ^ 
((7£^o(v))] with uncertainties in parentheses. An estimate of the noise in each data set is 
given as the width of the noise distribution, o j^^.^^{v); devices 2, 6 and 7 exhibit a 
correlation between enhanced o^^i^^iv) and larger o{v). Devices 1, 2, 5, and 7 (3, 4, and 
6) have a sheet density ~ 2 x 10^' cm~^ (3x10" cm~^) and mobility 
^-1.4x10^ crnlVs (6.5x10^ crn IVs). 
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